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Abstract. We present the result of photo-ionizing modelling of the three planetary nebulae (PNe) A 20, A 15 and 
MeWe 1-3. All three objects are roughly roundish, highly excited and have a high galactic z. The PNe displayed 
low densities in the shell, but relatively dense halos. A low metallicity and a relative high electron temperature 
were found. Comparisons with radio observations confirmed the obtained properties. The objects are very likely 
originating from thick disk stellar progenitors. The distances found investigating the PNe shells are somewhat 
lower than those derived spectroscopically with the central stars in the past. 
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\ 1. Introduction 

Although the interacting wind model for the formation 
of planetary nebulae (PNe) is accepted since a long time 
(Kwok et al. H978fl . the exact mechanism is not clearly un- 
derstood up to now. Several numerical simulations have 
been carried out especially during the last decade (e.g. 
Frank & Mellema l|1994afl . Dwarkadas & Balick 0998) 
or Schonberner & Steffen 1|2002|) . The last one carried 
out an one-dimensional hydrodynamic simulation based 
on the evolutionary tracks of Blocker ( 1995) and the simu- 
lations of the circumstellar shell around Asymptotic Giant 
■ Branch (AGB) star (Steffen et al. 11555)1 . 
The results of their simulation showed that a shock- 
bounded ionized main PN shell moving supersonically 
into the ABG-material compressing the inner parts of the 
matter into a dense but thin shell. The unaffected AGB- 
material becomes ionized as well and forms a rapidly ex- 
panding halo. Due to the drop of luminosity of the central 
stars of the PNe (CSPN) towards the white dwarf regime 
the outer part of the PN shell recombines quickly forming 
an second inner halo. The disadvantage of this simulation 
is that it just considers one dimension. Thus more complex 
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interactions of the winds like Rayleigh- Taylor instabilities, 
quickly building a clumpy environment as we see it in HST 
images, are not taken into account. This causes that the 
filling factor (e) was set in fact to unity. Comparing densi- 
ties obtain from forbidden line ratios and the one obtained 
from radio data or similar methods clearly show a filling 
factor far below unity (e.g. Mallik & Peimbert 1988] and 
Boffi & Stanghellini 1994). This changes the dramatically 
the optical thickness for ionizing UV photons and thus a 
recombining PN-shell is doubtful (Armsdorfer et al. 2002). 
The three PNe A 20, A 15 ( Abell IT555)> and MeWe 1-3 
(Melmer & Weinberger fIW0| were selected as they are 
roughly roundish objects. The central stars of the PNe 
(CSPNe) are well studied spectroscopically (Rauch et 
al. H555l McCarthy et al. 115571 Saurer et al. H557)l . This 
allows us to fix the input of radiation for the nebula. We 
modelled the physical properties of PNe-shells of three ob- 
jects and compare the results with the theory. 

2. Input Data 

2.1. Observations and Data Reduction 

Most data were obtained 1992 January 29 th and January 
30 th using the New Technology Telescope (NTT) with 
the EMMI spectrograph at ESO La Silla, Chile and were 
taken from the ESO data archive. THX31156 (red arm) or 
TK1024AF (blue arm) CCDs and grating #8 were used. 
The slit-width was adjusted to L"2. The spatial resolution 
was 0'.'44 per pixel at the red arm of the instrument (480 
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Fig. 1. The observations and the model for several lines of A 20. Left: Observation with NTT/EMMI (slit N-S); Right: 
Observation with ESO 3.6 m/EFOSC (slit E-W); especially in the [OIII] lines the radial symmetry of the object is 
well pronounced. Also the fragmentation in three parts, halo, PN shell and central cavity, is displayed well. The [Nil] 
lines are in an upper limit and thus confirm clearly the low metallicity found here (see Tab.[5J). 



nm to 680 nm) and 0'.'37 at the blue arm (400 nm to 505 
nm) . Additionally A 20 was observed by one of us (TR) in 
1996 February 10 th with the ESO 3.60 m telescope using 
the EFOSC1 spectrograph and the #26 TEK512M 1215- 
3 detector at a resolution of 0'.'61 per pixel (from 380 nm 
to 530 nm). At EFOSC observations narrow band [OIII] 
and H a images was taken too. 

The calibration was done using usual procedure in MIDAS 
and the calibration data from EMMI manual and from 
EFOSC1 manual. The standard star EG 21 (Hamuy et al. 
1992) was used for flux calibration. 

The target selection for this investigation followed two cri- 
teria: First, the PN should have a roughly round shape; 
second, the CSPN should be spectroscopically studied. 



Table 1. Basic data for the PNe sample: 



Name 




GPN 


Fl.4 GHz 


IRAS 




(Kimcswcngcr 2001 1 


[mjy] 




A 20 


G214 


96+07.81 


10.4 


07203+0151 


A 15 


G233 


53-16.31 


< 2.5 


06249-2520 


MeWe 1-3 


G308 


26+07.79 




13249-5426 



For all three objects IRAS counterparts have been 
found and radio-observations at 1.4 GHz have been found 
for two of them (Condon & Kaplan 1998 ). The objects are 
summarized in Tab. ^ 

2.2. Spectroscopy 

All three objects are characterized by an intensive Hell 
line (A = 468.3 nm) and a total depletion of [Nil] and [SII] 
lines. The [OIII] lines at A = 495.9 nm and A = 500.7 nm 
are much weaker than expected as well. The extinction was 
derived from the B aimer- decrement and to deredden the 



Table 2. Interstellar foreground extinction from the 
Balmer decrement derived here and from literature: 



Name 


E(B-V) 


E(B-V) 




here 


Tylenda et al. fTWZl 


A 15 


(TOO 


(F04 - (T64 


A 20 


(T05 


(F30 - m 35 


MeWe 1-3 


ff"27 





frames we used the interstellar extinction law by Savage 
k Mathis (fT5?9")> . The results are given in Tab.0 For A15 
and A20 our results from the Balmer-decrement of the 
high S/N spectra differ significantly from those derived 
from the survey spectra in Tylenda et al. (1992). After 
dereddening with our values from Tab. [21 the continuum 
of the CSPN fits perfectly to our results. Also the B — V 
of the CSPN from literature supports our results. The re- 
sulting (B — V)o is then always between -0™34 and -0™25. 
This is expected for the hot CSPNe. The relative errors 
of the line ratios were estimated by the variations of the 
ratios of Ha/H/3 and [OIII]5oo.7/[OIII]495.9 along the ra- 
dius. A conservative estimation gives an error of 10-15% 
for the relative line strength. Since the three [OIII] lines 
at A = 500.7 nm, A = 495.9 nm & A = 436.3 nm are de- 
tected a temperature determination was possible (Tab.[3J). 
The error of the ratio of [OIII] lines was assumed to be 50 
%, because the [Oiling 3 was very faint. Hence the accu- 
racy of T c was about T^™ K. Although the error of T c 
is relative high it turned out, that the electron tempera- 
ture is significantly higher than the 10 000 K, typical for 
PNe. This is a signature of the low metallicity (see below), 
causing the cooling by forbidden lines to fail. 
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Table 3. Electron temperatures derived from [OIII] lines 
and computed in CLOUDY: 



Name 


-^495. 9+^500. 7 
^436.3 


Te{K] 


(T e CL0UDY ) [K] 


A 15 
A 20 

MeWe 1-3 


54.7 
71.5 
62.2 


16 800™ 
14800™ 

15 800™ 


23 000 
20 000 
22 500 



2.3. Central stars 

The central stars of the objects have been studied very 
well in the past. Detailed spectroscopies were carried out 
by Rauch et al. ijlM^jl (A 20), McCarthy et al. l(TM7|> 
(A 15) and Saurer et al. ifTMTjl (MeWe 1-3). To determine 
distances Dcspn the method of Heber et al. 1|1984|) was 
used for A 20 and A 15 in the papers above. In the case 
of MeWe 1-3 no distance was published before. We thus 
did it using an Eddington flux H v given by Rauch ( 2004 ) 
following the method of Heber et al. (1984), who gave 
calibrations for Johnson and Stroemgren filters. 



Table 4. 



Name 


T e ff 


L 


log g 


Dcspn 




[kK] 


[Lo] 


[cm/s 2 ] 


[kpc] 


A 20 


119 ±22 


2100 


6.13 


4.1 ±50% 


A 15 


110 ± 10 


4238 


5.7 


4.5 ± 25% 


MeWe 1-3 


100 ± 10 


4580 


5.5 


6.9 ± 50% 



Unfortunately no reliable photometry is available for 
the CSPNe of A 20 and MeWe 1-3, thus the error is con- 
siderably higher. But also for A 15 an error of 25 % was 
assumed according Mendez et al. (1988 ). 

3. Photoionisation Model 

To model the PNe the CLOUDY code (Ferland lTM>j> was 
used. As input-parameters the radial density-profile, the 



chemical abundances and the stellar properties (within 
the ranges of error given above) were varied. The fill- 
ing factor was assumed to be 0.1, which is typical for 
evolved PNe (Mallik & Peimbert ifTMHjl and Boffi & 
Stanghellini lfT9"Mjl V The NLTE central star models of 
Rauch l|1997l I2003p were used to obtain a proper ioniz- 
ing spectrum. Armsdorfer et al. I|2002|) showed the need 
of real NLTE stellar atmospheres to describe especially 
the Hell lines properly. 

This results in the fraction of ionization of hydrogen 
and helium as a function of radius and a mean electron 
temperature. To obtain the surface-brightness, integration 
over the emissivity along the line of sight was carried out. 
To do this a spherical symmetry was adopted, what is in- 
dicated by round shape of the objects (see Soker (Q.997) 
and Tylenda et al. PTOljl for A 20 and A 15). The surface- 
brightness of various emission lines have been compared 
with the observation to achieve the best result. The av- 
erage electron temperature given by CLOUDY was com- 
pared with those derived before by the [OIII] lines. Due 
to the weak 436.3 nm line a spatial resolved study wasn't 
possible here. 

The objects A 15 and MeWe 1-3 displayed enhance- 
ments in the southern regions, especially in [OIII]. This 
might be due to interaction with the interstellar medium. 
Thus this regions where not taken into account for the 
model. 

3.1. A 20 (GPN G214.96+ 07.81) 

A 20 is the only one of these objects that is observed 
with both telescopes, NTT and ESO 3.60 m. The advan- 
tage of the doubled observation is that the slit direction 
was north-south at NTT and east-west at the ESO 3.60 m 
telescope and thus we were able to rule out systematic ef- 
fects in the data reduction and due to the assumed spher- 
ical symmetry. 

In Fig.nj&|2]the result of the modelling can be seen. The 
model was optimized with respect to the NTT spectra. 
In the east- west direction the Hell line is slightly stronger 
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than in north-south but both observed spectra agree fairly 
well with the model. Just in the very inner part the model 
displays an emissivity, especially of the [OIII] line, that 
is too high compared with the observations. This defi- 
ciency is, although uncertain, most likely originating in 
deviations from the spherical symmetry assumed in the 
CLOUDY code. The modelled [SII] and [Nil] lines are far 
below the limit of detection and thus confirmed by the 
observation. EFOSC-images do not display the regime, 
where the [Nil] and [SII] lines are located, but the line 
of [Nelll] at A = 386.9 nm and of [ArlV] at A = 471.1 nm. 
The [ArlV] line is very faint and the [Nelll] at the blue 
limit of the spectrograph, where the efficiency of the CCD 
and the transparency of the spectrograph have not good 
performances. Thus their error is higher than the one of 
the other lines. 
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Fig. 3. Density profile of A 20. The hydrogen- and the 
electron density are computed by CLOUDY. The heavy 
particle density is calculated from the hydrogen density 
and the chemical composition. 



In Fig. U the density profile of A 20 can be seen. The 
heavy particle density is mostly hydrogen and helium. 
Other elements contribute less than 2 • 10~ 4 . The high 
electron density, which exceeds the hydrogen density twice 
as much as the heavy particle density, indicates a mostly 
doubly ionized helium, what is confirmed in Fig. 

The chemical composition (Tab. 01 shows clearly an 
underabundance of heavy elements by a factor of 10. This 
makes it most likely to a member of the galactic thick disk 
population. The total depletion of [SII] and [Nil] lines can 
be explain just partial by the underabundance of these ele- 
ments. It is also caused by the fact that these elements are 
mainly in a higher ionization level. The overabundance of 
He is typical for PNe after the dredge-up processes during 
the RGB/AGB phase. 

The halo of the nebula is clearly highly ionized. It is 
not a recombination halo as predicted by ID hydro simu- 
lations. 
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Fig. 4. The fraction of doubly ionized helium as a function 
of radius. As can be seen helium is mostly doubly ionized, 
as it was expected considering the high Hell 468.6 nm 
intensity. 

Table 5. Chemical composition of A 20. The abundances 
of each element Z are given in log p^j and the comparison 
with the solar abundances (Grevesse et al. 2996) is given 
in dex. 





abundance 


abundance compared 
to solar abundance 


dominating 
ionization level 


He 


-0.825 


0.175 


He III 


C 


-4.45 


-1.00 


C IV 


N 


-5.00 


-0.97 


N IV 


O 


-4.10 


-0.97 


O IV 


Ne 


-4.824 


-0.90 


Ne IV 


Mg 


-5.523 


-1.10 


Mg III 


Si 


-5.523 


-1.07 


Si III 


S 


-5.824 


-1.15 


S IV 


Ar 


-5.80 


-0.32 


Ar IV 



3.2. A 15 (GPN G233.53- 16.31) 

In Fig. the surface brightness distributions for several 
lines of A 15 are displayed. There is hardly any decrease 
of the intensity towards the center of the nebula, contrary 
to A 20. This indicates that the inner part of the nebula 
is not empty, like in A 20, but somehow filled as can be 
seen in Fig. The intensity of [SII] and [Nil] is again, 
in conformity with the observations, below the threshold 
of detection. Also here the Hell line looks like the H/3. 
Another remarkable fact is, that the [OIII] is less promi- 
nent than the H/3 and the Hell line. 

In the inner part of the nebula the density relates to the 



radius as 



-0.3 



The main shell is thinner and has a lower 



density than the one of A 20 and followed by a relative 
extended and dense halo. The hydrogen in A 15 is also 
totaly ionized as well as most of the helium (Fig.EJ). 

The fraction of heavy elements is in this object even 
slightly lower than in A 20 and the electron temperature 
thus is higher as the cooling by the forbidden lines fails 
(see Tab. 01 what still agrees with the T e =23000K com- 
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Fig. 5. Comparison of the observation and the model of several lines of A 15. Ha and [OIII] (A = 500.7 nm) are not 
taken into account since they are related to H/3 and [OIII] (A = 495.9 nm) respectively. The [Nil] and the [SII] was 
not observable. The innermost 'hole' in the data is faked due to the CSPN subtraction/overlay. 
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Fig. 6. Density profile and fraction of doubly ionized helium of A15. Description like in Fig. |3& Fig. 0] 



puted by CLOUDY). Also the dominating ionization lev- 
els are slightly higher than the one in A 20, especially the 
one of C, S, and Si. The typical He enrichment is shown 
by A 15 too. The overabundance of nitrogen with respect 
to the other heavy elements give us not an hint for a mas- 
sive progenitor, but originates in the upper limit of the 
observed line. 

3.3. MeWe 1-3 (GPN G308.26+07.79) 

MeWe 1-3 shows, even if it is not so pronounced like the 
one in A 20, an inner central cavity. Also in respect to 
other conditions it looks like A 20, as there are the similar- 
ity of the Hell and the H/3 line, the relative line strength 
of [OIII] and the lack of any [SII] and [Nil] lines. Only 
there is hardly any halo visible in MeWe 1-3. 
Since the emissivity distribution of MeWe 1-3 looks like 
the one of A 20 their density profile is similar - an empty 
inner part, and a rapid increase of the density to a plateau, 
here at 100 H cm" 3 . Just the narrowness of the halo of 
only 0.013 pc is different. In this case the whole nebular 



Table 6. Chemical composition of A 15 (same as in 
Tab.0. 





abundance 


abundance compared 
to solar abundance 


dominating 
ionization level 


He 


-0.82 


0.18 


He III 


C 


-4.95 


-1.50 


C V 


N 


-5.03 


-1.0 


N IV 


O 


-4.60 


-1.47 


O IV 


Ne 


-5.420 


-1.50 


Ne IV 


Mg 


-5.920 


-1.50 


Mg III 


Si 


-5.950 


-1.5 


Si V 


S 


-6.17 


-1.5 


S V 


Ar 


-6.98 


-1.5 


Ar IV 



hydrogen and helium is mostly bare of all electrons too. 
Also the abundances of the elements are very similar to the 
one of A 20. But the average ionization stage is closer to 
the one of A 15. Also the electron temperature (15 800 K 
observation, 22 500 K model) is in between the two other 
nebulae, and that despite the lower T g of its central star. 
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Fig. 7. Comparison of the observation and the model of several lines of MeWe 1-3 (description see Fig. [5J . 
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Fig. 8. Density profile (left) and fraction of doubly ionized helium of MeWe 1-3. 



The larger luminosity of the CSPN of MeWe 1-3 and the luminosity and the distance, we varied them within the 
smaller radius of the nebula itself is able explain this ef- ranges of error. The obtained data are given in Tab. [H] 
feet. 



Table 7. Chemical composition of PN MeWe 1-3 (same 
as in Tab. 01. 





abundance 


abundance compared 
to solar abundance 


dominating 
ionization level 


He 


-0.85 


0.15 


He III 


C 


-4.45 


-1.0 


C IV 


N 


-5.03 


-1.0 


N IV 


O 


-4.13 


-1.0 


O IV 


Ne 


-5.03 


-1.11 


Ne IV 


Mg 


-5.42 


-1.0 


Mg III 


Si 


-5.45 


-1.0 


Si V 


S 


-5.67 


-1.0 


S V 


Ar 


-6.48 


-1.0 


Ar IV 



Table 8. Obtained basic data of the PNe modeled. 



Name 


Ten [kK] 


L [L e ] 


radius [pc] 


dist. [kpc] 


A 20 


140 


3000 


0.40 


2.35 


A 15 


120 


4200 


0.32 


4.01 


MeWe 1-3 


110 


4500 


0.21 


3.95 



It seems that the temperature and thus the luminosity 
of the CSPN of A 20 are underestimated. The uncertain 
distances of A 20 and MeWe 1-3 are always smaller in the 
models compared to the one given in the literature, but 
still within the error-deviation. 



4. Discussion and Conclusion 



The metallicities of all three PNe investigated are lower 
Since the errors of the properties of the objects, given by a factor of 10 to 30 compared to the solar abundances, 
in the literature, are relatively large, especially for the only helium is little bit more present than in the Sun. 
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Soker ( 2002 ) already reports an under-abundance of met- 
als in round PNe. 

Looking at the galactic coordinates and the distances of 
the planetary nebulae one can easily verify that the height 
above the galactic plane (z) of these nebulae is compara- 
tively large (Fig. |5J . 

The nebulae investigated are located at a z of 320 pc, 
-1127 pc and 535 pc for A 20, A 15 and MeWe 1-3 re- 
spectively. To compare these values with other one the 
PNe catalogued by Acker et al. (1992) were taken as a 
representative sample for all PNe. Half of these planetary 
nebulae are located below a z of 162 pc. Their average z 
amounts to 288 pc. The large z of the three objects inves- 
tigated indicates that they are so called thick disk objects 
and therefore their progenitor was an underabundant ob- 
ject as well. 
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Fig. 9. Here the galactic height z versus the radius from 
the galactic center is displayed. The radius of the sun was 
assumed to be 8 kpc. The + shows the PNe catalogued 
by Acker et al. (fT59l)l and the the PNe A 20, A 15 and 
MeWe 1-3. 



The high electron temperature of the nebulae is also 
linked to the low metallicity (Stasinska |1 978| ) . As Frank 
& Mellema l|1994bj> showed one important cooling mech- 
anisms in ionized gas is the radiation emitted by colli- 
sionally excited lines. The most important role plays the 
[Nil] lines in low excited nebulae and the [OIII] lines in 
high excited one. Since no [Nil]- and less intensive [OIII] 
lines, compared to common PNe, are detected a higher 
temperature is feasible. 

As mentioned above a filling factor (e) of 0.1 was 
adopted. Since radio data are available for A 20 and A 15 
(Condon & Kaplan I1998f) it was possible to check this 
assumption. Adopting a optically thin wind in the radio 
regime and therefore \ v oc v~ 01 allows us to compute the 
radio flux at 5 GHz from the 1.4 GHz data and thus the 
total ionized mass according Cahn et al. I|1992|l . 



M ion . = ^2-266 x 10- 21 • D 5 • 6> 3 • F 



5 GHz 



where D is the distance in pc, the angular radius in arc- 
seconds and F 5 qhz the radio flux at 5 GHz in Jansky. We 
compared this mass with the one of our model, computed 

by 

r+00 

M = 47T/1 / n hp (r)r 2 dr 
Jo 

ji is the mean atomic weight and nh p (r) the heavy par- 
ticle density as function of the radius. This results in 
e. Although the error of the obtained filling factor is 
quiet high it is clearly smaller than 0.3. Using the sam- 
ple and the formula of Mallik & Peimbert (1988) we ob- 
tain a value even below 0.05. The formula of Boffi & 
Stanghellini (1994) give a value of 0.096. But especially 
the sample shows no PN with a radius above 0.2 pc and 
an filling factor above 0.1. Tests with CLOUDY show that 
the nebulae are always optically thin for the UV radiation 
up to e = 0.5. 

Further calculations showed that even stars with lumi- 
nosities down to 160 Lq are able to ionize most of the 
hydrogen and thus a recombination of the shell seems to 
be impossible. This is in clear contradiction to the pre- 
dictions of the ID hydro simulations (e.g. Schonberner & 
Steffen EOOIJ ■ 

While the luminosities of the CSPN found from the 
investigation of the PN shells here are about the same as 
those found with the NLTE models of the CSPN spec- 
troscopy. Our temperatures are - still within the error 
bars - slightly higher than those from the CSPN spec- 
tra. As shown by Lechner & Kimeswenger 1|2004|) and 
Emprechtinger et al. ( 2004 ) these parameters are very sen- 
sitive to the line ratios of the highly ionized species. We 
thus suggest this effect to be real. The distances Dqspn 
have to be significantly lower than those found by other 
investigations in the past. 
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